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1. Introduction 

When thinking about quadratic and homogeneous Hamiltonian systems on the 
dual of a Lie algebra, the first example that comes to our mind is the system 
describing the rotations of a free rigid body around its center of mass. This system 
was derived by Euler in 1758 (see e.g. [5]) then generalized by Poincare (see [14]), 
and then again later by Arnold (see e.g. [2]) starting from the original so(3) Lie 
algebra, to a general Lie algebra. The mathematical literature contains a huge 
amount of writings concerning Euler's equations, from their original form to the 
most general forms (see e.g. [5], [12], [1], [3], [4], [6], [7], [8], [15], [10], [13], [9], [11]). 

Recently, in [17] it is proved that under suitable non-degeneracy conditions, a 
quadratic and homogeneous (quadratic) Hamiltonian system defined on the dual of 
the Lie algebra o(K) of real K - skew-symmetric matrices, where K is an arbitrary 
3x3 real symmetric matrix, is afimcly equivalent to Euler's equations (Euler's 
equations with three linear controls) describing the free rigid body dynamics. More 
exactly, in both cases above, the non-degeneracy condition requires the existence 
of a positive definite linear combination between the quadratic form generated 
by K and respectively the quadratic form that generates the Hamiltonian of the 
system. In [17] it is introduced a general family of quadratic Hamiltonian systems 
that generalizes the quadratic Hamiltonian systems defined on the dual of the Lie 
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algebra o(K), and consequently gives also a generalization of Euler's equations 
describing the free rigid body dynamics (with linear controls). 

The purpose of this paper is to construct explicitly the normal form of an 
interesting subset of quadratic Hamiltonian system from this general family, that 
generalize Euler's equations (with linear controls). 

For details on Poisson geometry and Hamiltonian dynamics, see, e.g. [1], [2], 
[13], [6], [7], [8], [16]. 

2. Normal form of a large class of three dimensional quadratic 
Hamiltonian systems 

In this section we give an explicit method to compute the normal form of a general 
quadratic Hamiltonian system on a natural extension of (o(K))* , where K G 
Sym(3) is a given 3x3 real symmetric matrix. 

Recall that the Lie algebra o(K) is defined as o{K) := {A G g[(3; R) : A T K + 
KA = 03}. Note that in the case when K is nonsingular, o(K) is the Lie algebra 
of the Lie group 0(K) := {A G GL(3;R) : A T KA = K} of K-orthogonal 3x3 
matrices. 

Let us now recall a Lie algebra isomorphism between o{K) and R 3 . For details 
regarding this isomorphism, see e.g. [6], [13]. 

Proposition 2.1. The Lie algebras (o(K), +, -r, [•, •]) and (R 3 ,+,-r, Xjf) are iso- 
morphic, where [•,•] is the commutator of matrices, and respectively uXk%) := 
K(u x v), for any u,v G R 3 . 

Note that for K nonsingular, the above isomorphism between S G o(K) and 
s£l 3 , can be defined by using the equation Su ~ s x Ku, iieR 3 . 

An immediate consequence of the this proposition is that (o(K))* = (R 3 )* = 
R 3 , viewed as a dual of a Lie algebra, it has a natural Poisson structure, namely 
the "minus" Lie- Poisson structure, which in this case proves to be generated by 
the Poisson bracket: 

{/, 9 } A -:=-VC x .(V/xV 3 ), 
for any f,g G C°°(IR 3 ,R), where the smooth function C K G C°°(R 3 ,R) is given by 

C K (u) := ^u T Ku. 

For more details regarding this bracket and the associated Lie-Poisson dynamics, 
see e.g. [6]. 

It is not hard to see that the center of the Poisson algebra C°° (R 3 , R) is 
generated by the Casimir invariant function Ck G C°°(R 3 ,R), Ck(u) = —u T Ku. 

Let us now generalize the above setting by introducing the Poisson manifold 
(R 3 , {•, •}nri t )), where the Poisson bracket {•, -}(km) is defined by: 

U,9}(KM) :=-VC m -(V/xVj), 
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for any f,g £ C°°(R 3 ,R), and the smooth function C {K k) £ C°°(R 3 ,R) is given 

by 

C(K,k)(u) ■= -u T Ku + u T k. 

Using the above definition of the Poisson bracket, it follows that the smooth func- 
tion C(x,k) is a Casimir invariant function. 

Consequently, a quadratic Hamiltonian system on (R 3 , {■, -}(^ k )), is gener- 
ated by a smooth function i?(A,a) G C°°(R 3 ,R), given by 

H( A ,a.){u) ■= ^u T Au + u T &, 

where A £ 5ym(3) is an arbitrary real symmetric matrix, and a £ R 3 . 

Hence, the associated Hamiltonian system is governed by the following dif- 
ferential equation: 

u = (Ku + k) x (Au + a), u £ R 3 . (2.1) 

The main result in [17] states that if there exists a, (3 £ R such that aA + (3K 
is positive (or negative) definite, then the system (2.1) is affinely equivalent to 
Euler's equations of the free rigid body dynamics with linear controls, namely: 

u = u x (Du + d), u £ R 3 , (2.2) 

where D = diag(Ai, A2, A3) is a real diagonal 3x3 matrix, and d £ R 3 . In the case 
when d = 0, the equations (2.2) are exactly the Euler's equations of the free rigid 
body dynamics. Note that in coordinates, the system (2.2) become: 

ii = (A 3 - A 2 ).T 2 a;3 + d 3 x 2 ~ d 2 x 3 , 
&2 = (Ai - A 3 ).xia;3 - d 3 xi + dtx 3 , 
X3 = (A2 - \\)xix 2 + d 2 xi - dix 2 , 

where (di,d 2 ,d 3 ) are the coordinates of d. 

The purpose of this paper is to give a normal form of the equations (2.1) in 
the case when the matrices A and K commutes. Let us now state the main result 
of this paper. 

Theorem 2.2. If the matrices A and K commutes, then the system (2.1) is orthog- 
onally equivalent to the dynamical system: 

v = (D K v + k) x (D A v + k), v £ R 3 , (2.3) 

where the real diagonal 3x3 matrices Da, Dk are given by D A = R T AR, D K = 
R T KR, where R £ 0{ld) = 0(3, R) is a 3 x 3 orthogonal matrix, and k := 
det(R)R T k, a := det(i?)i? T a. 

Proof. As a consequence of the commutativity of the real symmetric matrices A 
and K (i.e., [.A, if] = AK — KA — 0), there exists an orthogonal matrix R £ 
O(Id) = 0(3, R) which diagonalize simultaneously both the matrices A and K. 
More exactly, there exists two 3x3 real diagonal matrices Da and respectively 
Dk such that Da = R T AR, and respectively Dk = R T KR. 
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Let us denote R T := (R 1 ) T . Recall that for any matrix R £ O(Id) = 
0(3, R) we have that (det(i?)) 2 = 1. 
Next, we will show that: 

u(t) = det(R)Rv(t), 
where t M- u(i) is a solution of the dynamical system (2.1): 

u = (Ku + k) x (Au + a) , 
and respectively t <— >• v(t), is a solution of the dynamical system 

v = (Dkv + k) x (Dav + a). 
To prove this assertion, note first that: 

(Ku + k) x (Au + a) = (K det(R)Rv + k) x (Adct(i?)i?w + a) 

= [det(R)KRv + (det(i?)) 2 k] x [det(R)ARv + (dct(i?)) 2 a] 

= (dct(R)f[(KRv + dct(i?)k) x (ARv + dct(i?)a)] 

= (KRv + dct(i?)k) x (ARv + det(i?)a) 

= (RD K v + det(JJ)k) x (RD A v + det(R)a) 

= (RD K v + det(R)RR T k) x (RD A v + det(R)RR T &) 

= [R(D K v + dct(i?)i? T k)] x [R(D A v + dct(i?)i? T a)] 

= det(R)R~ T [(D K v + k)x (D A v + a)] 

= det(R)R[(D K v + k) x (D A v + a)]. 

Hence, 

u = (Ku + k) x (Au + a) dct(R)Rv = dct(R)R[(D K v + k) x (D A v + a)] 

v = (D K v + k) x (D A v + a). 

□ 

Remark 2.3. Using coordinates, the system (2.3) becomes: 

±i = (K 2 A 3 - K 3 A 2 )x 2 x 3 + (K 2 a 3 - k 3 A 2 )x 2 + (k 2 A 3 - K 3 a 2 )x 3 + k 2 a 3 - k 3 a 2 , 

1 2 = (K 3 Ai - KiA 3 )xix 3 + (k 3 Ai - A'ia 3 )iri + (K 3 ai - /ciA 3 )x 3 + k 3 ai - fcia 3 , 

1 3 = (KiA 2 - K 2 Ai)xix 2 + (Kia 2 - k 2 Ai)x\ + (k\A 2 - K 2 a 1 )x 2 + k x a 2 - k 2 a\, 

where D A = diag(Ai, A 2 , A 3 ), D K = diag(A'i, K 2 , K 3 ), a = (ai,a 2 ,a 3 ), and k = 
(fci,fc 2 ,fc 3 ). 
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